DO-branes with non-zero angular momentum 

George K. Savvidy 



National Research Center Demokritos, Ag. Paraskevi, GR- 1 5 3 10 Athens, Hellenic Republic 
Abstract. 

In my talk I shall consider the mechanism of self-expansion of a system of N DO-branes into high- 
dimensional non-commutative world-volume investigated by Harmark and Savvidy in [ 1 ] . Here D2- 
brane is formed due to the internal angular momentum of DO-brane system. The idea is that attractive 
force of tension should be cancelled by the centrifugal motion preventing a D-brane system from 
collapse to a lower-dimensional one. I shall also present a new extended solution where a total of 
9 space dimensions is used to embed a DO-brane system. In the last section, by performing linear 
analysis, the stability of the system is demonstrated. 



INTRODUCTION 

In the last years there has been increasing interest in dimensionally reduced supersym- 
metric Yang-Mills theories [2, 3, 4]. One of the reasons is that the reduction of ten- 
dimensional theory to p + 1 dimension is relevant for the description of Dp-branes, p- 
dimensional extended objects carrying Ramond-Ramond (RR) charges in type II super- 
string theories [5]. In the extreme reduction to zero dimensions it is believed to describe 
DO-branes, fundamental pointlike objects in type IIA superstring theory [6]. In certain 
energy regimes the dynamics of N such particles can be described by the supersymmetric 
quantum mechanics ofNxN Hermitian matrices obtained from dimensional reduction 
of 3\£ = 1, D = 10 super- Yang-Mills theory down to + 1 dimensions 1 . 

It is believed that supersymmetric quantum mechanics of many DO-branes in type IIA 
superstring theory is equivalent to a partonic description of light-front M-theory [1 1], a 
more fundamental underlying M-theory [8, 9, 10]. The existence of matrix formulation 
of M-theory [11], the BRSS-conjecture, crucially relies on the existence within type-TYA 
string theory of a tower of massive BPS particles electrically charged with respect to 
the RR 1-form. These particles originally described as black holes in IIA supergravity 
were identified with DO-branes and can be interpreted as Kaluza-Klein particles of 
eleven-dimensional M-theory compactified on a circle. The existence of the M-theoretic 
Kaluza-Klein tower of states is equivalent to the statement that supersymmetric Yang- 
Mills quantum mechanics has exactly one bound state for each N [12, 13, 14]. 

Remarkable aspect of matrix theory is that not only classical gravitational interactions 
can be produced in the large N-limit, but also the appearance of the superstring extended 
objects in terms of pointlike fundamental degrees of freedom. One example of such 



It was also proposed that + 0- dimensional matrix model should give a Poincare invariant non- 
perturbative definition of IIB superstring theory, the so-called IKKT model [7]. 



phenomena is the observation that a IIA superstring and a system of N DO-branes can 
be blown-up to a D2-brane by placing it in background field [15, 16]. 

In my talk I shall consider another mechanism of self-expansion of a system of N DO- 
branes into high dimensional non-commutative world-volume investigated by Harmark 
and Savvidy in [1], here D2-brane is formed due to the internal angular momentum of 
DO-brane system. The idea is that attractive force of tension should be cancelled by the 
centrifugal motion preventing a D-brane system from collapse to a lower-dimensional 
one 2 . I shall also present a new extended solution where a total of 9 space dimensions 
is used to embed a DO-brane system. In the last section, by performing linear analysis, 
the stability of the system is demonstrated. 

In the D-brane formulation the dynamics of N DO-branes can be described by the 
supersymmetric quantum mechanics of N x N Hermitian matrices obtained from di- 
mensional reduction of fA£ = 1, D = 10 super- Yang-Mills theory to 0+1 dimensions 
[6, 16, 20, 21] (the quantum mechanical model was originally studied in [4, 22, 23]). 
The effective action of N DO-branes is the non-abelian SU(N) Yang-Mills action plus 
the Chern-Simons action 

S YM = -To(2nlj) 2 JdtTr(^F, v F^ , (1) 

where is the non-abelian SU (N) field strength in the adjoint representation and 
To = (gJs)^ 1 is the DO-brane mass. To write this action in terms of coordinate matrices 
X 1 , one has to use the dictionary [1] 

with i,7= 1 , 2, . . . , 9 and in Aq = gauge we have 

Sym = T q J ^TrQx^-i^^[r,^][xU^ . (3) 
The equations of motion are 

and should be taken together with the Gauss constraint 

[X',X ! ]=0. (5) 

The Chern-Simons action derived in [16, 20, 21] for the coupling of N DO-branes to bulk 
RRCW andd 3 ) fields is 



Scs = T J dtTr(c + QX i +^ [c$[X i ,X J ]+C$[X i ,Xi]X k 



(6) 



2 Similar phenomena appear in the cases of rotating branes on spheres [17, 18], "giant gravitons", and in 
the cases when angular momentum is generated by crossed electric and magnetic BI fields [19]. 



and describes the interaction of ./V DO-branes with slowly varying background fields of 
Type IIA supergravity. Even though the DO-brane world- volume is only one-dimensional 
a multiple DO-brane system can couple to a brane charges of higher dimension. Myers 
[16] considers the system of N DO-branes in a constant external 4-form RR field strength 
F4 = dC^ and has found a stable solution, where the DO-branes are polarized and 
arranged into a static spherical configuration. A lower-dimensional object under the 
influence of higher-form RR fields may nucleate or be 'blown-up' into D2-brane. The 
background field imposes an external force that prevents the collapse of the D-brane to 
a lower-dimensional one. 

Another way to get self-support against collapse is to allow DO-brane system to carry 
mechanical angular momentum [ 1 ] . This new kind of rotating solution of the system of N 
DO-branes was constructed by Troels and Konstantin in [1] and in the subsequent articles 
[24, 25] it was demonstrated that this solution describes a stable D2-brane configuration. 
Below I shall concentrate mostly on this solution and on it generalizations presented in 
[25]. The basic idea in their construction is that the attractive force of tension should be 
cancelled by the centrifugal repulsion force. The earlier work where membrane solution 
appeared with non-zero angular momentum was [26]. 



D2-BRANE FROM MULTIPLE ROTATING DO-BRANES 

I shall briefly review the spherical D2-brane configuration of type IIA string theory since 
a new solution for a rotating system of N DO-branes presented below uses the essential 
elements of this construction. The solution is equivalent to the spherical membrane 
solution of M(atrix) theory [1 1, 29]. 

With the aim to construct a membrane with an S 2 geometry we shall embed the S 2 in 
a three-dimensional space spanned by the 123 directions and consider the ansatz 

Xi(t)= J— Un(t), 1 = 1,2,3 (7) 

VN 2 - 1 

where the N x /V matrices Li , L2, L3 are the generators of the N dimensional irreducible 
representation of SU (2), with algebra 

[L;,Ly] = iEijkLk • (8) 

and with the quadratic Casimir Yd=i = > so mat Tr(L 2 ) = Af ^ 2 ~ 1 \ For vanish- 
ing background fields the Hamiltonian is H = ^ \ £? =1 ff + ^ (r\r\ + r\r\ + r\r\ 

where we have introduced a convenient parameter a = y ^i- This gives the equa- 
tions of motion 

f x = -a 2 (rl + rj)r h r 2 = -a 2 (r 2 + rj)r 2 , r 3 = -0C 2 (r 2 + rj)r 3 . (9) 

This system is otherwise known as + 1 dimensional classical SU(2) YM mechanics 
[22, 23]. Let us for simplicity take all radii to be equal to each other: r\ = r 2 = r$ = r. 



With this we have from (7) the physical radius of the membrane R 2 = Xf +X% +X 2 = 
I r 2 , where I is the N x N identity matrix. The last formula shows that the /V D0- 
branes are constrained to lie on an S 2 sphere of radius r. The equations of motion (9) in 
this case reduce to the equation r = —2a 2 r 3 and the solution is r(t) = Rq sn(Ro<xt + §) 
oscillating between Rq and — Rq 3 . One can trust this solution if \r(t) | <C l s VN, \r(t) | <C 
1, \r(t) \ <C l~ . Since we also require \r(t)\ ^> l s we must have N >> 1. Thus it is 
necessary in order to have a large amount of DO-branes to build a macroscopic spherical 
membrane. 

In order to construct the rotating ellipsoidal membrane, viewed as a non-commutative 
collection of moving DO-branes we shall take previous configuration of the non- 
commutative fuzzy sphere in the 135 directions, and set it to rotate in the transverse 
space along three different axis, i.e. in the 12, 34 and 56 planes. We thus use a total of 6 
space dimensions to embed our DO-brane system. The corresponding ansatz is [1] 

x 1 (t) = -^=Un(t) , x 2 (t) = -^=Ur 2 (t) , 

X 3 (t) = -^=Ur 3 (t) , m = 7 ^=jUr 4 (t) , 

X 5 (t) = - / ^=L 3 r 5 (t) , X 6 (t) = -^=L 3 r 6 (t) . (10) 

where the N xN matrices Li , L2, L3 are the generators of the TV-dimensional irreducible 
representation of SU(2). In this ansatz the matrix structure is such that the coordinate 
matrices are proportional to the SU (2) generators in pairs and the Gauss constraint (5) 
is identically satisfied. It is also the only finite-dimensional subalgebra of the group of 
diffeomorphisms of S 2 , the Sdiff(5 2 ) [28]. That is why the SU (2) ansatz is in some sense 
unique: it is the only type of solution that carries over to the supermembrane without 
modification [30]. 

Substituting the ansatz into (3) gives the Hamiltonian H = ^ ^ £f =1 r] + ^ {r\ + 

r i) ( r 3 + r \) +( r i + r i) ( r 5 + r l) + ( r 3 + r 1) ( r 5 + r i) ) and tne corresponding equa- 
tions of motion 

h = -a 2 (rl + rl + rj + rl)n , h = -a 2 (rf + r\ + rj + r\) r 2 , 
r 3 = -a 2 (r 2 + rj + r 2 5 + r\) r 3 , r 4 = -a 2 (rf + rl + rj + r\) r 4 , 
r 5 = -a 2 (r 2 + rl + r 2 ,+rl)r 5 , r 6 = -a 2 (r 2 + r 2 2 + r 2 3 +r|) r 6 . (11) 

The special solution of these equations, describing a rotating ellipsoidal membrane with 
three distinct principal radii R\,R 2 and R 3 is [1] 

n(t) = i?icos(cflif + (|>i) , r 2 (t) =i?isin(0Di? + (|)i) , 



3 Clearly the spherical membrane will be classically point-like at the nodes of the sinus. Thus the 
membrane solution will break down after a finite amount of time, since the classical solution may not 
be valid at substringy distances, and possibly decay into a Schwarzschild black hole [29]. 



n(t) = R 2 cos((£> 2 t + § 2 ) , r 4 (t) = R 2 sin((£> 2 t + § 2 ) , 

r 5 (t) = R 3 cos((D 3 t + § 3 ) , r 6 (t) = R 3 sm(& 3 t + 4> 3 ) • (12) 



This particular functional form of the solution ensures that the highly non-linear equa- 
tions for any of the components r,- are reduced to a harmonic oscillator. The solution (12) 
keeps r\ + r\ = R\ , r\ + r\ = R\ and r\ + r\= R 3 fixed, which allows us to say that the 
object described by (12) rotates in six spatial dimensions as a whole without changing 
its basic shape. At any point in time one can always choose a coordinate system in which 
the object spans only three space dimensions. 

Using the equations of motion (11), the three angular velocities are determined 

by the radii, and do not necessarily have to coincide: (Oi = a^R\-\-R 2 , a>2 = 

OLtJr 2 + R 2 , (0 3 = &\Jr\ + R 2 ■ This dependence of the angular frequency on 

the radii is such that the repulsive force of rotation has to be balanced with the attractive 
force of tension in order for (12) to be a solution. Thus the radii R\,R 2 and R 3 parame- 
terize (12) along with the three phases §i, to produce altogether a six-parameter family 
of solutions. In order to exhibit the properties of the solution (12) one shall evaluate 
the energy, the components of angular momentum and to find out nonzero D-brane cur- 
rents. The non-zero components of M t j are M\ 2 = —M 2 \ , M34 = — M43 and = —M^ 
and correspond to rotations in the 12, 34 and 56 planes respectively. Their values fit 
with the interpretation of the solution as N DO-branes rotating as an ellipsoidal mem- 
brane in that they are time-independent due to conservation law and proportional to N: 
Mn = \NTQto X R\ , M 3A = \NT Q ® 2 R 2 , M 56 = \NT Q (H 3 R 2 , thus 

M 2 = ±(NT ) 2 (w 2 l Ri + & 2 2 R 4 2 + w 2 3 R 4 3 ) , E = ^ (co 2 * 2 + co 2 ^ + co 2 ^) . (13) 

Let us now compute the coupling of the solution (12) to the RR potential. The 
interaction with the C^-field is governed by the action (6). We denote the corresponding 
current by J l i k 

/J " = 2^ Tr ( [X[, ' XJ] ^ ] ) ' (14) 

and we have to impose anti-symmetrization with respect to ijk indices. For our solution 
the non-zero components of Q are: <2i35 = ^NaRiR 2 R 3 cos(X>itcos(X> 2 tcos(X> 3 t, together 
with 2 2 46 ,<2l46 , 2i36 , Q235 , 0236 , 2i45 , £?245- From this it is easy to obtain the 
corresponding 7's by differentiation with respect to time. Thus the Chern-Simons action 
(6) shows that the coupling of this system to F0123 is non-vanishing and that the spherical 
membrane solution has a D2-brane dipole moment. The higher currents are equal to zero 
and our system does not carry D4-D8-charges. 

In addition I shall present "breathing" brane solutions. For that it is convenient to 
introduce polar coordinates (p cos(|), p sintf)) so that the Hamiltonian // = takes the 

form: H = \ Yfi=\ [pf + P 2( P] + 2~[pfp2 + P2P3 + P3Pi]-The conservation integrals are: 
pf cf>i = Mi, p2<[>2 = M 2 , p 2 § 3 = M3, where M,- = ^ and the Hamiltonian takes the 



form: H = ±Z? =1 



■2 , M} 



+ ^[Pip2 + P2P3 + P3P1] • The equations of motion are 



Pi = -p ; (p 2 - pj) +Mf/pf , 1= 1,2,3 



(15) 



and our previous solution (12) is p, = 7?,- = Const, i= 1, 2, 3 and <j>? = to? = Mf/Rj = 
R 2 +l + R} +2 - In the special case when all coordinates are equal to each other pi = 
P2 = p3 = p(0 and can depend on time, the "breathing" brane solution, then H = 
I [P 2 + + P 4 ] an( l corresponding equation can be integrated. The new solution is 

elliptic function p = p(t) [24] 4 . 

In order to increase the number of parameters of the rotating N DO-brane system, I 
shall take previous configuration and set it to rotate in the transverse spaces along three 
different axis, i.e. in the 123, 456 and 789 planes. Thus I shall use a total of 9 space 
dimensions to embed DO-brane system. The corresponding ansatz is 



X 1 (t) = 

x 4 (t) = 

X 7 (t) = 



2 

2 



Un(t) , 




L2T4(0 , 






Xg(t) 



2 
2 



Vn 2 



Ur 2 (t) ,X 3 (t) = 
L 2 r 5 (0 ,X 6 (t) = 
L 3 r 8 (0 ,X 9 (0 = 



2 
2 



Lir 3 (?) , 
L2 r6(0 > 

L 3 r 9 (/016) 



The coordinate matrices are again proportional to the SU (2) generators and 
the Gauss constraint (5) is identically satisfied. Substituting the last ansatz into 

(3) gives the Hamiltonian: H = ^f- (\Y? i = l r 2 + ^ {r\ + r\ + rf) (rj + r\ + 
rf) + (rj + r 2 + r 2 )(r 2 + r 2 + r 2 ) + (r 2 + r 2 + r 2 )(r 2 + r 2 + r 2 ) ) and the equations 



n = -a 2 {r\ + r 2 5 +rl + r 2 1 + rl + rj) n, i = 1 , 2, 3 

fj = -a 2 (rj + r 2 2 + r 2 3 + r^ + rl + r\) rj , j = 4,5,6 

== ™2/„2 , „2 , „2 , „2 , „2 , „2 

^5 " 



-CC z (rf + ti + rj + ri + rf + ri)r k , fc = 



7,8,9. 



(17) 



The special solution of these equations, which ensures that the highly non-linear equa- 
tions for any of the components are reduced to a harmonic oscillator, is: 

n(t) =/?icos((Dif + <|>i) , r 2 (t) =Ri sin(0i) • sin(cQi? + ,r 3 (t) =i?icos(0i) - sin(coif + <j>i) , 
r 4 (t) = R 2 cos((0 2 t + § 2 ) , rsCO =#2sin(e 2 ) -sm((Q 2 t + § 2 ) ,r 6 (t) = R 2 cos(Q 2 ) ■ sm((d 2 t + <h) , 
H(t) = R 3 cos(& 3 t + § 3 ) , r 8 (0 = # 3 sin(e 3 )-sin((0 3 / 1 + 4> 3 ) ,r 9 (t) = R 3 cos(B 3 ) ■ sm(& 3 t + 0$) 

The object described by (18) rotates in nine spatial dimensions as a whole with- 
out changing its basic shape. The radii R\, R 2 and R 3 parameterize (18) along with 



4 A somewhat similar ansatz was proposed in [31], where some of the features of the solution were 
foreseen. 



the six phases 0, and fy, to produce altogether a nme-parameter family of solu- 
tions. The energy and the angular momentum are the same (13). The interaction 
with the C^ 5 ) -field is governed by the Chern-Simons action derived in [16, 20, 21] 
j^fdtSTr(c$ kl [X\Xj][X k ,X l ] +clf klm [X i ,XJ][X k ,X l ]X^ . We denote the corre- 
sponding current by J l i klm 

jiju™ = —}—STr([X^,xj}X k [X l ,X m ^]) . (19) 

One can be convinced that all components of J are equal to zero. Thus the Chern-Simons 
action shows that the coupling of this system to F^j^im is vanishing and that the extended 
solution does not carry D4-brane charge. 



STABILITY ANALYSIS WITHIN SU(2) AND FULL SU(N) GROUP 

The purpose of this section is to present a complete stability analysis of the fluctuations 
in the neighborhood of the rotating DO-brane solution of [1]. Initially in [1] were 
analyzed perturbations that do not modify the original SU (2) ansatz. In [25, 24] this 
analysis was extended to the case when perturbations are in the full SU (N) algebra 
directions. In the full SU(N) case there are exactly A^ 2 + 12 zero-modes, of which 
N 2 — 1 are the consequence of the global color rotation symmetry of the solution, and 
6 are associated with global space rotations. All other modes are completely stable and 
execute harmonic oscillations around the original trajectory. 

Let me present the stability analysis of the SU (2) (I = 1) perturbation of the system 
[1] and then move to more general SU(N) (I = 2,3,4,....) perturbations [25, 24]. The 
equations of variation which follow from equations in polar coordinates (15) are: 

8pi = -4(J?2 + J^)8pi -2KiJ? 2 5p2-2KiJ? 3 5p3 (20) 

502 = -2i? 2 #i8pi -4(R 2 + R 2 3 )8p 2 -2R 2 R 3 8p3 (21) 

5p 3 = -2/? 3 *i5pi -2fl 3J K 2 8p2-4( J R 2 + J R 2 ! )Sp3, (22) 
and have only positive modes [24] 

Q. 2 = 4( R 2 +R 2 +R 2 ), (23) 

a\ 3 = 2(R 2 + R 2 +R 2 ) ± ^2((R 2 -R 2 ) 2 + (R 2 -R 2 ) 2 + (R 2 -R 2 ) 2 ). (24) 

To consider perturbations in all directions of underlying SU (N) group we should repre- 
sent SU (N) generators Y l m as higher order monomials in the N x N matrix generators 
L/, i = 1,2,3 of the SU(2) group 

Y L= Y< C m (ii,...,i z ) L »i ' ' ' • ( 25 ) 

h,...,ii 

The total number of generators is then £^7/(2/ + 1) = N 2 — 1 as it should be for 
SU(N). The general explicit construction of the Y l m is due to Schwinger (see also 
[32, 33, 27, 34]). 



Let us represent the rotating DO-brane solution in a more convenient form: X'(t) = 
LjR cos(otf), X l (t) = X i+3 (t) = L t R sin(atf), z = 1,2,3. In what follows we will set 
R=l, with (O equal to u) 2 = 2i? 2 = 2. In the basis provided by the spherical operators 
Y l m (25) we have [28, 30], 



T Y' 



mY 



for /= 1,...,N-1 



L±Xl = v/(/ T m)(/±m+l)Y^ ±1 . 



(26) 



We will not use the explicit form of these matrices, as the defining relations (26) is 
all that is needed. The properties under Hermitian conjugation can be summed up as 
V Jm = ( — l) m Y/,_ ra . I shall first consider the perturbations which are parallel in space 
to one of the directions of our solution, that is in the directions 123456 and then in the 
directions 789. Let us decompose the fluctuation fields in the basis defined by Y l m , where 
/ runs from 1 to /V — 1 



5X' = £ yU' 



m=—l 



5X = £ y^C 

m=—l 



i= 1,2,3. 



(27) 



We do not explicitly show an / index on the T|, ^. We should also impose 
(-l) m £_ m and r)^ = (-l) ra r|_ m for all m = —I,. ...I . The variational equa- 
tions of motion are 



-5X< = [8X' , [X', X 1 ] ] + [XJ , [SKi, X 1 ] ] + [XJ , [X', 8X'] ] + 
[8# , [V, X'] ] + [tJ , [5X>, X'] ] + [V , [V, 5X'] ] . 

The constraint equation looks like 



£ [SX',X'1 + [x',8X'j + [5X',X'1 + [x<,5X' 



= 0. 



(28) 



(29) 



Using the commutation relations (26) we get for the constraint ^,LJ ;m (cos((0?)4m + 
(Dsm((Ot)^ l m + sm((Dt)i] l m — (Qcos((tit)r\ l m ) =0 , where L^ m are now the 5/7(2) generators 
in the (2/ + 1) x (21 + 1) representation. In the co-moving coordinates 



4 = cos((Qf)& + sm((QfK, and = -sin((Qf)^, + cos(ajf)Tii, (30) 



the constraint looks simpler, m U nm (u l m — 2(Ov l m ) = . (31) 

The variational equation of motion (28) after substituting the fields (27) is 

g + 1(1 + 1) ^ = cos((Qf ) (L^L^ + le^L*,,) (cos((Qf ) ^ + sin((Qf ) t^) . (32) 

The decoupling of the modes with different / is seen to be a direct consequence of 
(26), and more fundamentally, of the pure SU (2) structure of the original background 



solution (12). The equation for T] is gotten by exchanging cosines for sines and £ for T). 
In the co-moving coordinates (30) the equation becomes a linear system with constant 
coefficients: 

u l n + (l(l+\)-2)u l n -2av l n = (Li n ,Ll, m + iE jik L k nm )ui, (33) 
v;+(/(/ + l)-2)v^ + 2co4 = 0. (34) 

Thus we shall analyze the system of equations (31) (33), (34). The rhs of (33) is a matrix 
acting on a 3(2/ + 1) component vector and the eigenvalues A of this block matrix are 
given in the table, together with their multiplicity [25]. Choose a fixed frequency ansatz 
u l n (t) — e l0x u l n , v l n (t) = e lQt v l n . The second equation (34) can be solved as 

-2^2 iQ 2 

V " = /(/ + l)-2-Q^ - (35) 

and then substituted back into the first equation (33). In the basis in which the ma- 
trix on the rhs is diagonalized, it can be replaced with its respective eigenvalue 

A, resulting in an algebraic equation for the Q.: (l (I + 1) — 2 — Q. 2 ) —8£l 2 = 
A (l (I + 1) — 2 — Q. 2 ) . Finally, this quadratic equation can be solved, Q 2 2 = 

— jA + 1(1 + 1) + 2± j a/ A 2 — 16A + 32/(/ + 1) and the corresponding modes are 
given in the table 



A 


Q 2 


a 2 


multiplicity 


1(1 + 1) -2 





l 2 +l + 6 


2/ + 1 


21 


/ 2 -3/ + 2 


l 2 + 3l + 2 


2/ + 3 


-(21 + 2) 


I 2 -I 


l 2 + 5l + 6 


2/-1 



(36) 



Note that the number of zero modes changes from 9 for the case / = 1 and 12 for 
/ = 2, to 21 + 1 for arbitrary / > 2. Thus the total number of zero modes is the sum 
9 + 12 + (21 + 1) = N 2 + 12. Now I shall consider perturbations that are in the 
directions 789, if we had oriented the original solution along 123456. The perturbations 
§X k = Yirn^m^m f° r £ = 7, 8,9 satisfy the simple harmonic equation + /(/ + 
1)^ = . This clearly has only positive frequencies and is therefore stable. For / = 1 
all the 9 modes have the same frequency as the original solution, corresponding to 
infinitesimal global rotations of the system into the 789 hyperplane. The counting goes 
as follows, there are 9x2=18 first order degrees of freedom here, which coincides with 
the dimensionality of the grassmanian manifold of embeddings of a 6-hyperplane into 

^' e ' 50(6)x5Q(3) • F rom these results it follows that, zero-modes notwithstanding, all 
the frequencies in the system are positive, and arbitrary small perturbation will remain 
bounded for all times 5 . 



5 The same problem was considered also in the paper [35]. However the authors of [35] initially arrived 
at the Mathiew equation instead of the equations (33), (34), (31) and therefore to the opposite conclusion, 
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